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Abstract. The Jaynes-Cummings (JC) model is a milestone in the the- 
ory of coherent interaction between a two-level system and a single 
bosonic field mode. This tutorial aims to give a complete description 
of the model, analyzing the Hamiltonian of the system, its eigenvalues 
and eigestates, in order to characterize the dynamics of system and sub- 
systems. The Rabi oscillations, together with the collapse and revival 
effects, are distinguishing features of the JC model and are important 
for applications in Quantum Information theory. The framework of cav- 
ity quantum electrodynamics (cQED) is chosen and two fundamental 
experiments on the coherent interaction between Rydberg atoms and 
a single cavity field mode are described. 



1 Introduction 

The JC model was originally proposed in 1963 by Edwin Jaynes and Fred Cummings 
PQ in order to study the relationship between the quantum theory of radiation and the 
semi-classical theory in describing the phenomenon of spontaneous emission. In the 
semi-classical theory of atom- field interaction, the field is treated as a definite function 
of time rather than as an operator, whereas the atom is quantized. The semi-classical 
theory can explain many phenomena that are observed in modern optics, for example 
the existence of Rabi cycles in atomic excitation probabilities for radiation fields 
with sharply defined energy (narrow bandwidth). The JC model aims to find how 
quantization of the radiation field affects the predictions for the evolution of the state 
of a two-level system, in comparison with semi-classical theory of matter-radiation 
interaction. What is a purely quantum effect, described only by the JC model and 
not by the semi-classical theory, is the revival of the atomic population inversion after 
its collapse, providing a direct evidence of the discreteness of field states (photons). 

The JC model is still an important and actual topic in quantum physics, since 
it touches different branches of the most recent and outstanding research activities, 
and big efforts have been carried out to extend this basic theoretical topic [3J. The 
natural framework for the implementation and experimental testing of the JC model 
is cQED, where fundamental experiments has been performed using microwave and 
optical cavities 4 5 6 2. The JC model has been also used in the framework of trapped 
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ions, where the stable elctronic states of the ions interact with their vibrational modes 
(phonons) [7]. In the context of solid state systems, semiconductor quantum dots are 
placed inside photonic crystal, micropillar or microdisk resonators |5]. Nowadays the 
advanced technology in circuit-QED allows to build a two-level system from quantum 
circuits interacting with transmission line microwave resonators [3]. 

In order to more precisely describe the interaction between an atom and a laser 
field, the model is generalized in different ways. Some of the generalizations com- 
prehend different initial conditions jfOj . dissipation and damping in the model 
multilevel atoms and multiple atoms [12 , and multi-mode description of the field [T3"] . 
Other extensions of the JC model include a driving pump laser acting on one or more 
two-level atoms, which allows the simultaneous action of rotating and counter-rotating 
coupling terms in the Hamiltonian, including dissipative effects due to Markovian en- 
vironments [14115116117] . 

The JC model serves as a tight connection between Quantum Optics and Quantum 
Information Theory, since it naturally provides the dynamics of interacting qubits, 
the basic amount of quantum information encoded in every two-level system. Several 
JC-based systems have been employed in order to perform quantum computation 
algorithms, quantum teleportation and quantum cryptography protocols |18j . 

The present tutorial is structured in a pedagogic form, in order to provide all the 
mathematical passages necessary to the comprehension of the basic JC model and its 
figures of merit. The physical model is introduced in Sec. 1, then the uncoupled and 
coupled systems are analyzed respectively in Sec. 2 and Sec. 3. The diagonalization 
is performed in Sec. 4, bringing to the complete diagram of the energy levels of the 
system. Sec. 5 is devoted to the description of the dyanamics, with a particular stress 
on the resonant case, where the transition frequency of the two-level atom is equal 
to that of the cavity field mode. A formalization of the JC model is performed in 
Sec. 6, by introducing the density operator of the system and of the subsystems. 
Vacuum Rabi oscillations and phenomena of collapses and revivals are the objects of, 
respectively, Sec. 7 and Sec. 8. The final part of the tutorial is devoted to illustrate 
two model experiments in cQED, where the quantum nature of matter-radiation 
interaction shows itself in a clear and powerful way. 



2 The model 



The fully quantized JC model concerns a two-level system interacting with one ra- 
diation mode. Under the dipole approximation and the rotating wave approximation 
(RWA), this kind of system is described by the following Hamiltonian: 

H=H + Uaf = -^ r 3 + M flt a + \) + %(ar+ + oV") (1) 

where Haf = h~g{ar + -\-a)r~) is the interaction term and Ho is the free Hamiltonian. 
The atomic transition frequency is oja, the frequency associated to the field mode is 
u), while the coupling constant g quantifies the strength of the interaction between 
the atom and the electromagnetic field. The atomic excited state is denoted as |e), 
with energy E e = Huja/2 and the ground state as \g), with energy E g = —huiA/^- The 
atomic raising operator is r + — \e)(g\, the atomic lowering operator is r~ = \g)(e\ 
and r3 = |e)(e| — \g)(g\ is called inversion operator. The atomic transition energy 
associated the inversion operator r% is Hlua — E e — E g . The field mode creation 
operator is a\ whereas the field mode annihilation operator is a. 
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The operators introduced above satisfy the commutation relations 

[a,a f ] = I 

[r + ,r-] = r 3 (2) 
[r 3 ,r±] = ±2i* 

providing the following conditions for the Hamiltonian operators: 

— at resonance uj = lo a 

\Rq,W.af\ = and also [%o,%] = 
meaning that "Ho and % have common eigenvalues; 

— the total number of excitation quanta is conserved: [N, Ti\ = 0, where N = a^a + 
r + r~ is a constant of motion. 

2.1 Rotating wave approximation (RWA) 

It must be noted that the original and general Hamiltonian describing the interaction 
of a two-level system with a quantized field mode, is different from Eq. (TTJ) and 
includes the so called anti-rotating terms, r + aJ and r~a. In fact, under the dipole 
approximation, the two-level atom can be described by a dipole moment al coupled 
to an electromagnetic field E oc (a + a^) and the interaction reads 

Hi=-d-E. (3) 

The atomic dipole moment can be expressed in terms of the atomic basis states and 
only the off-diagonal terms are non-zero since, by parity considerations, (e|c?|e) = 
(<? \d\ g) = and one can write 

d = d\g){e\+d*\e){g\=d(r- +r + ) (4) 

where (e\d\g) = d and it has been assumed, without loss of generality, that d is real. 
Hence, the interaction between the two-level atom and the quantized field mode is 

Ui = hg{r+ +r~)(a t + a) (5) 

where the coupling constant g depends on the strength of the interaction between the 
atomic dipole and the field mode. The Hamiltonian describing the whole system is 
then % = Ho + T-Li- This is the general way to describe the quantized matter-radiation 
interaction, which goes under the name of Rabi model. 

In order to recover the Hamiltonian (JTJ) describing the fully quantized JC model, 
the RWA is employed to discard the anti-rotating terms, r + a* and r~a. This fact is 
appreciated if the system dynamics is described in the interaction picture, considering 
the transformation Itf'(i))/ = Ut\^(t)) with XJ t = exp(j:Hot) and where the state 
vector \^(t)) is a solution of the Schrodinger equation 



ih-Mt)) = u\*(t)). 



(6) 
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A similar expression can be derived for the state vector \W(t))i in the interaction 
picture 

m^\nt))i=Ki{tmt))i (7) 

with 

Hj{t) = U}UiU t = fig (r + ah i ^ +UA)t + rae-' l{ul+UJA ^ + r+ae" 1 ^-'"^* + raV^"^') 

(8) 

Close to the resonance condition, uj ~ cj^, the terms r + a^ and ra are multiplied by 
the factors qM^+^ap which are fast rotating in time. The RWA allows to neglect these 
terms in describing the dynamics of the system and to recover the JC Hamiltonian 
(P), since 

Hi{t) ~ fi. 9 (r + ae-' l{u - UA)t + raV^"^*) (9) 

which gives in the Schrodinger picture the typical JC interacting term Haf- 

A more mathematical approach for the RWA is to consider the Dyson series expansion 

of the time evolution operator in the interaction picture 

Ui(t)=I-±J 'Hi(jt 1 )dti + (-£\ J dhJ X Hi(t 1 )H I (t 2 )dt 2 + ... (10) 
which rules the time evolution of the state vector in the interaction picture 

\*(t))i = u I (W(o))i- (ii) 

Considering, for instance, the first order in the series expansion, the integration pro- 
vides terms as ±(w + uja)^ 1 connected to the operators r + a) and ra, whereas terms 
like ±(w — uja)" 1 are linked to r + a and ra) . The last terms are dominant with respect 
to the first ones under the resonance condition uj ~ uja, together with the assump- 
tion that <C 1. With these conditions, the major contributions to the time 
evolution operator come from the resonant terms r + a and ra) . 

The JC model is valid only if the RWA can be employed to describe an effective 
coherent interaction between a two-level system and a field mode. This can be done 
whenever the coupling constant g is much smaller than the single frequencies uj and 
uja- Recently, new regimes has been explored both theoretically |19j and experimen- 
tally [2UJ where the coupling constant g is comparable with the atomic and field 
frequencies, for which the RWA is no longer applicable. 



3 Uncoupled system 



Without any interaction between the atom and the cavity field mode (g — 0), the 
Hamiltonian of the uncoupled system is 



n = n = h 



r 3 + uj I a} a + 



1 



(12) 



One can distinguish two main situations, the resonant and the non resonant case: 



At resonance (uj — uja) Eg. IT21 becomes 

H = huj + a) a + = huj(r + r- + a f a) = HujN (13) 
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Finding eigenstates and eigenvalues of "Ho is a quite easy task. Referring to the 
previous observation about the conservation of the number of quanta of the sys- 
tem, one can consider the system states \e,n) = |e) \n) (atom in the excited 
state and n photons) and \g 7 n + 1) = \g) (g) \n + 1) (atom in the ground state and 
n + 1 photons), where |n) (with n = 0,1,...) denotes the Fock basis states of a 
bosonic field mode. These states are eigenstates of the number operator N with 
the same eigenvalue n + 1 and, hence, of the system Hamiltonian Hq with the 
eigenvalue hw(n 4- 1). The only exception is the ground state of the whole system 
\g, 0), for which H \g, 0) = holds. 



In general, for a non-resonant condition, the Hamiltonian of the uncoupled system 
can be written as 

'r 3 + 1\ hS 



(14) 



where the detuning parameter S = uja — w has been introduced. Now, as depicted 
in Fig.[TJ the two system eigenstates \e, n) and \g, n + 1) have different eigenvalues 



H \e,n) = E e , n \e,n) = h 
H \g 1 n+ 1) = E g<n+1 \g,n+ 1) = H 



(n+l)uj + - 



(n + l)u 



e,n) 
5' 



\g,n + l) 



(15a) 
(15b) 



giving rise to n orthogonal Hilbert subspaces ^(n) of dimension 2x2. 



4 Coupled system 

Introducing the coupling term in the Hamiltonian of the system, Haf — hg(ar + + 
aV~), the previous separated systems (atom and field mode) become a unique en- 
tangled system. One can notice immediately that Haf couples only those states that 
belong to a single subspace (or bidimensional manifold) S'(n). In fact: 

— the process of absorption is mathematically described by the matrix element 



(e, n\HAF\g, n + 1) = hg\/n + 1, 

which corresponds to the transition \g, n + 1) — > |e, n); 
— the emission process is given by 



(g,n+ l\H AF \e,n) = hgyjn+ 1, 

which corresponds to the transition |e, n) \g, n + 1). 

Hence, it is sufficient to consider only the bidimensional manifold, eigenspace of Hq 
at a given number of photons n: 

S(n):{\e,n),\g,n + 1)} (n = 0, 1, 2, . . .). (16) 

Within each of these subspaces, the Hamiltonian of the system is given by the matrix 

(g,n + l\H \g 1 n + 1) (g,n+ l\HAF\e,n)\ 

(e,n\'HAF\g,n+ 1) (e,n\H \e,n) I 

(17) 



(n) n 
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Fig. 1. Scheme of the naked levels -E e ,n and E g: „+i of the uncoupled system described by 
Ho. Each couple of energy levels corresponds to a couple of vector eigenstates, which form a 
bidimensional Hilbert subspace S(n). The ground state \g, 0) is independent from the others. 



which can be rewritten as 



(H)„ = Hu)(n+l) 



1 
1 



(18) 



where Q n = 2gy/n + 1 is called Rabi frequency for n photons. 



5 Diagonalization 



In this section the eigenvalues and eigenvectors of Eq. ([18]) are calculated. 
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5.1 Eigenvalues 

By imposing the eigenvalue equation 
/(n + l)w + f 
V % (n + lK 

one obtains the two energy eigenvalues: 



2 




(n + l)w± 



ft* 



(19) 



where 7?.„ = v/ J? 2 + <5 2 is called generalized Rabi frequency from which, at resonance 
conditions 6 = 0, one gets the Rabi frequency lZ n = fl n . At resonance the degeneracy 
of the two energy states of the system is lifted and the energy separation is given by 
(AE) n = HQ n . In this sense the Rabi frequency plays the role of a Bohr frequency for 
the whole system and the energy separation of the two levels is higher as the number 
of photons increases, (AE) n oc \fn (see Pig. 



5 > 



E e 



E, 



n A 

hs 

n+l V — 



.(re 4- l)hw- 



5 = 



5 = (resonance) 



(n + 1)Huj 



9 = 



± E 



Fig. 2. Scheme of the dressed levels E±„ of the Hamiltonian T-L of the coupled system, 
related to that of the uncoupled system. At resonance 5 = the degeneracy is lifted by the 
coupling and the separation between the two energy levels is (AE) n = hf2„. 



5.2 Eigenstates 

The eigenstates of the system corresponding to the eigenvalues (fT9|) are also called 
dressed states, as the "naked" atomic states are now "dressed" by means of the in- 
teraction with photons. By imposing the eigenvalue equation in matrix form (|18[) one 
obtains 

5 fl n \ ( u \ . ( u 



and after some passages, using the normalization condition u 2 + v 2 = 1, the eigen- 
vectors are found to be: 



1 / a 



y/(K n - 6) 2 + Ql \K n - 6 
u\ If lZ n — S 



vi ~ y/(n n - s) 2 + ni v -a 



(21a) 
(21b) 
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The two eigenvectors can be rewritten in a more compact form, in terms of basis 
eigenvectors, as: 



10 



COS I 



" 10 



COS f 



sin f 



(22a) 
(22b) 



and tan 29 n = - ^ 



where sin 9 n = , , cos 9 n = 

Going back from the eigenvectors (matrix form) to the state vectors (Dirac notation) 
within the manifold £{n), one has: 



H\±n) =E± jn \±n) =h 



(n + l)u) ± 



Tin 



±n) 



where the dressed states in fucntion of the naked states are 



- n) = sin & n 
n) = cos 9 n \e, n 



|e, n) + cos9 n \g,n + 1) 
— sin 9 n \g, n + 1) 



(23) 



(24) 



which are general entangled states of the global system atom-field. 
The interaction Hamiltonian Haf produces a unitary transformation on the basis 
states {|e,n), \g,n + 1)}, which is simply a rotation of an angle 9 n in S(n). In partic- 
ular, at resonance condition one obtains: 



(25) 



which is a maximally entangled state, with 9 n = j. From the relation tan2(9„ = — 
one gets the states (j25|) as S -> 0, while, if |<5| — > co the mixing angle 9 n — > and the 
uncoupled (or naked) states are recovered. 

Within each subspace S"(n) and varying the detuning parameter 5, the corresponding 
levels' chart is depicted in Fig. [3] 




Fig. 3. Dressed eigenvalues in function of S = uja — IjJ. The coupling removes the degeneracy 
at S = 0, whereas the naked eigenvalues cross each other. When S 3> 1 the dressed eigenvalues 
tend to the naked ones, as the atomic and field subsystems were uncoupled. 
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5.2.1 Naked levels 

The eigenvalues of Hq are linear in <5: 



(n + l)w+- 
(n + l)w- - 



and they cross at S = showing a degeneracy at resonance E e ^ n = E gjl+ i = (ri+l)hui. 



5.2.2 Dressed levels 

As obtained before, the eigenvalues of % are E± n H [(n + 1)uj ± , with lZ n = 

At resonance (S = 0) the energy levels do not cross, but they are separated by the 
effect of the interaction by the quantity (AE) n = hQ n = 2hgy / n + 1. 
On the other hand, quite off resonance \S\ ^> fi n , the generalized Rabi frequency is 
lZ n ~ \5\ and the dressed states are asymptotically equal to the naked states. This 
phisically means that the eigenstates of the total Hamiltonian H coincide with those 
of the uncoupled Hamiltonian T-Lq or, equivalently, it means that g — > 0. 



6 Dynamics 

In this section the dynamics of an initial pure state under the action of the Hamilto- 
nian (TTJ) is considered. In the Schrodinger picture the evolved state is given by 

\V(t)) =e-* K *|#(0)) (26) 

where the initial state vector 1^(0)) can be expanded on the basis of the dressed states 
{| +n), | — n);n = 0, 1, . . .} as 

oo 

l*(°)> = E (j + n)(+nmo)) + I - n)(-n\V(0))) 

r (27) 

= ^(c +n (0)| + n)+c_ n (0)|-n)) 
with c± n (0) = (±ra|<P(0)). Inserting Eq. ((37J) in Eq. one obtains: 

OO 

!*(*)> = E (c-* B+ " t C+ n(0)| + n) + e-* B -*c-n(0)| - n)) (28) 
n=0 

where the eigenvalues E± n are given by Eq. (TT9"]) , 

Within a given eigenspace S'(n), the vector state at time t is 

|ff(t))„ = e - i ("+ 1 ^ t (e-^^c+^O)! + n> + e i ^" t c_„(0)| - n)) (29) 

which can be written in the interaction picture \&{t)) n ->• |^(*))« = e i(n+1)wt |ip r (t)) n , 
eliminating the free oscillation at the mean frequency (n + l)w, as 

(i) | + n) + c_„ (i) | - n) (30) 
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where the evolved coefficients are defined as c± n (t) = e =F1 T 2l 'c± n (0). 
Projecting Eq. (|3U)) on the basis states {{— n\, (+n\} one obtains the complex coeffi- 
cients c± n (t) in the matrix form 

fc_ n (f)\ _ U*** o \/c_„(oA 

^(tjj-^ o e-*>*y W„(0)J (31) 

In order to discuss about the physics of the dynamics, Eq. (I3T1) must be expressed 
on the naked basis {|e, n), \g,n + 1)}. The complex coefficients c± n (t), c e ^ n {t) and 
Cg,n+i(t) are related by the unitary transformation 

/cos^-sinO) 

where cos 9 n and sin (9„ have been defined previously. Equation pip becomes 

( C e , n (t) \_rrff \ f C e ,„(0) \ 

{c g , n+1 (t)J ~ U [ o e-^t) U [c g , n+1 (0)J 
which brings to the explicit relation 

/ Ce , n (t) \ = -^ sin (¥*) -i|: sin (¥*) \ / Ce ,„(o) \ 

^.^-iWj ^ _i|^ sin (^t) C os(^)+i^-sm(^)J U.«+i(°)y 

(33) 

that describes the time evolution of the complex coefficients of the naked eigenvec- 
tors. Physically one has the complete and general description of the dynamics of the 
whole system, in terms of the atomic energy levels and the number of cavity photons, 
including detuning, coupling strength and initial conditions for pure states. 



6.1 The resonant case 

At the resonance condition S — 0, the generalized Rabi frequency lZ n coincides with 
the Rabi frequency Q n — 1g\Jn + 1. In this case Eq. (|3"3"]l assumes the following simple 
form 

Ce,n(t) \ = ( ^(f^) - isiQ (^)\ / C e ,„(0) 

c g , n+ i(t)J ^_ isin (^^ cos(^i) J V c ^»+i(°)/ 

As an example one can consider the initial pure state |<^(0))„ = |e, n), that is, an 
excited atom is injected in a cavity with the field mode in a Fock state with n 
photons. The complex coefficients at t — arc c e ,„(0) = 1 and c ff ,„+i(0) = 0. After a 
certain interaction time t > 0, the evolved coefficients are 



(34) 



(35) 



.„(£) = cos(^i) 
c 5 , n +i(t) = -isin (^t) 

from which the evolved state of the system can be written as 

m))n = oos(^t) \e,n)-i S m(^t) \g,n+l). (36) 
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Fig. 4. Plot of Rabi oscillations for the probability amplitudes P e ,„ and P gi „+i in function 
of the dimensionless time f2 n t. The initial state of the system is |^(0)) = |e, n). 

From the relations one can calculate the probability amplitudes P e , n (t) = |c ej „(i)| 2 
and P g n +i{t) — \c g n +i{t)\ 2 corresponding, respectively, to the states {|e, n), \g, n + 

i n 

P e , n (t) = COS 2 (^t) = \ [1 + COB(«„t)] 

(37) 

P fll „+i(t) = sin 2 (^f) = i [1 - cos(tf„*)] . 

The levels occupation probabilities of the JC system oscillate at the Rabi frequency 
Q n = 2g\Jn + 1, frequency at which the two-level atom and the field mode exchange 
a single photon, process known as Rabi oscillations or Rabi flopping. fl n is, actually, 
the transition frequency for the dressed states, which are the energy levels to be 
considered in the JC model with interaction between the atom and the field mode. 
As plotted in Fig. [4] and using Eqs. (l35|) . at times f2 n t = kn (with k = 0, 1, 2 . . .) the 
state of the system, if k is even, is the basis eigenstate ±|e, n), or =Fi| (7, n + 1) if k is 
odd. At intermediate times fi n t = (2k + l)ir/2 the system is in an entangled state, 
as for instance -4= (|e, n) — i|<?, n + 1)) when Q n t = tt/2. 



7 Formalization 

In this section a more mathematical treatment of the time evolution of a generic 
state of the system is introduced, together with the analysis of the atomic and field 
subsystems. 



7.1 Evolution operator 

The time evolution operator of the JC model, in interaction picture, is given by 

n=0 



1 Trivially, if the initial state of the system is 1^(0)) = \g,n + 1), the probabilities (|37p 
become P e , n (t) = sin 2 ( ■%*■£) and P 9jIl +i(t) = cos 2 f^rf-i). 
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where the interaction Hamiltonian is Haf = hg(ar + + a'r ). Developing Eq. (|38l) 
one can divide the sum into even and odd terms like 

' (ar+ +a)r-) 2n = (ad) n \e)(e\ + (ala) n \g)(g\ 
(ar+ +a)r-) 2n+1 = {aa<) n a\e){g\ + a* {aa}) n \g){e\ 

resulting in a time evolution operator of the form 

(-i 5 i) 2 " 



Ujc(t) = Y,{- 



n=0 



+ 



(2n)! 



(Vaa^) 2n \e)(e\ + (Vafa) 2n \g)(g\ 



2n+l 



(2n + l)! 



i jJ= a\e)(g\ + =4=^- 



l<7)(e| 



that, rewritten in a more compact way, it appears as 



U JC {t) = cos{gWN + l)\e)(e\+cos(gtVN)\g)(g\ + 

t sin(gty^r+T)) 



1 ,__ a\e)(g\ - ia> 



VN + T 



VN + T 



\9)(e\ 



(39) 



(40) 



where the number of photons operator is N = a 1 a. Written in a matrix form the time 
evolution operator assumes the form 



(41) 



where u gg (t) = cos(gtVN), u ge (t) = -ia f sm( ^ =±^, u eg (t) = -i 81 "^^ — a and 
tt ee (t) = cos(gty/N + 1). 



7.1.1 Pure state 

A generic initial pure state for the atom-field system can be written on the naked 
basis {|e, n), \g,n)} (which includes also the ground state |g,0)) in the form 



mO))=J2(c e , n (0)\e,n)+c g , n (0)\g,n)) 



(42) 



with Cen(O) = c e (0)c„(0) and c gn (0) = c 9 (0)c„(0). The evolved state \&(t)) 
UjcWm is 



n=0 

c g (0) 



\#(t)) = f>(0){c e (0) cos^t) |e,n>-isin(^) \g,n + l) 

tj |e,n-l) } 



+ 



— 1 \ I , / f^n — 1 \ i 

cos I — — t \g, n) - ism — — t \e, n - 1) 



(43) 



As an example one could take as the initial state of the system |^(0)) = |e, n), 
verifying that the evolved state (j43]) is exactly that of Eq. (|36p . 
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7.1.2 Mixed state 

A generic initial mixed state of the system is described by the density operator p(0) = 
p F (0) <E> Pa{0), where the field and atomic subsystems can be considered factorized 
at t = Q. The time evolution of the system is 

P(t) = Uj C (t)p(0)ul c (t). (44) 

As an example on could take as initial state of the system the density operator 
p(0) = p F (0) <8> |e) (ej, where the atom is in the excited state and the cavity mode in 
a generic mixed state, and the evolved state of the system is 

p(t) = u ee (t)p F (0)ul e (t)\e}(e\ + u ge (t)p F (Q)u\ e {t)\g){g\+ 

+ + ( 45 ) 

+ u ee (t) PF (0)«t e ( t ) \ e ){g\ + u ge (t)p F (0)ul (t)\g) (e\. 



7.2 Atomic time evolution 

In order to consider only the atomic subsytem dynamics, one must apply a partial 
trace of the system density operator over the field degrees of freedom. Starting, as 
usual, from the simple case of p(0) = p F (0) ® | e) (e | and using Eq. (l4"5j) the density 
operator describing the atomic time evolution is 

p A (t) = Tr F [p(t)] = Tv F [u ee (t)p F (0)ulM\e)(e\+Tr F [u ge (t)p F (0)ul e (ma)(9\ 
+ Tt f [u ee (t) PF (O)^ e (f)] \ e ){g\ + Tr F [u ge (t)p F {Q)u\ e (t)]\g) (e \ . 

(46) 

For an experimental measurement the observables are the atomic probabilities, which 
can be expressed as 

> B (t) = (e\ PA (t)\e) = Tr F [ul(t)u ee (t)p F (0)] = (u\ e u ee ){t) 

(47) 

P g (t) = (g\p A (t)\g) = TY F [ul e (t)u ge (t) PF (0)} = {u\ e u ge )(t) 

where the cyclic property of the trace operation^ has been used. As can be easily 
veryfied, the sum of the probabilities must give one, 

P e (t)+P g (t) = Tr F [{4 e (t)u ee {t) + u\ e {t)u ge {t))p F {Q)] = Tr F [I F p F (0)] = 1 

A particular example can be considered, that of the field mode prepared in a diagonal 
form p F (0) = J2 m Pm\m){'m\i with the normalization condition ^2 m P m = 1 and 
where P m is the probability corresponding to the Fock state |m). This is the case 
of the field in a Fock state \n) with P m = S m ^ n , or in a thermal state with P m = 

( i+{N) ) ' w here (N) is the average number of thermal photons. In this case, 

the off diagonal terms in Eq. (|46p are null and the density operator for the atomic 
subsystem is a statistical mixture of the atomic states |e) and \g): 

p A (t)=P e (t)\e)(e\+P g (t)\g)(g\ (48) 

where the probabilities P e (t) and P g {t) are calculated in Eqs. (|4T|) . 



The cyclic property of the trace states: Tr[ABC] = Tr[CAB] = Tr[BCA] 
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7.3 Field mode time evolution 

The dynamics of the field mode may occur in two ways, depending on whether the 
state of the atom exiting the cavity is measured or not. In general the two types of 
measurement are labelled as "selective" and "non-selective" measurements which are, 
respectively, a non-unitary and a unitary evolution. 

If, at a particular time t, the atom is measured outside the cavity it has passed 
through, the field mode inside the cavity, which is entangled with the atom due 
to their previous interaction, undergoes a non-unitary evolution. The field state is 
projected into a normalized state, quantum correlated with the atomic state measured 
at that particular time. As before, one can consider as the initial state of the system 
p(0) = p_f(0)(8> |e)(e| and, depending on the measured atomic state |e) or \g), the field 
state is 



PF(t) 



f o (e) (t) - ""(t)PF(Q)«I (*) _ u ee (t) PF (o)ut (t) measured in le^ 

P F W ~ Tr F [u ce (t) PF (0)uUt)] ~ Pe(4) at ° m meaSUreCl 111 \ e ) 

p f(t) = "r (t) /T ( °i"rw vi = "- (t)p ;S° K ° (t) atom measured in | 5 ) 



On the other hand, if the atom is not measured or, equivalently, the outcome con- 
cerning the atomic state is ignored (non-selective measurement), the field state will 
be a statistical mixture of the two possible outcomes 

PF {t) = P e {t) P { ;\t) + P g {t)pf{t) = u ee (t)p F (0)ul(t) + u ge (t)p F (0)ul e (t) 
= F$(t)p F (0) 

where J 7 ^ \(t) is the "super-operator" of Jaynes-Cummings for the reduced dynamics 
of the field alone, with the atom initially in the excited state. This is also the way 
of writing a quantum map, that is an operation that transform a density operator in 
another density operator, leading to a physically acceptable description of the system 
evolution [T8:. Usually it is written as a Kraus sum 



P 



c{ P ) = m, p mi with K M * = 1 ( 51 ) 



where M M are called Kraus operators. As regards the system above, the Kraus op- 
erators are {u ee (t), u ge (t)}, while F e = u\ e (t)u ee (t) and F g = u ge (t)u ge (t) are called 
POVM (Positive Operator Valued Measure), which are generalized measurement op- 
erators that include the class of projectors Pj with the property pj = Pj. POVMs 
must satisfy only the property of summing up to one, that is F e + F g = I, which 
ensures to have positive probabilities as P e - g = Tv[pF e - g ]. 



8 Vacuum Rabi oscillations 

Retrieving the previous analysis on the Rabi oscillations for resonant condition S = 0, 
the particular case with the initial condition |^(0)) = |e,0) can be considered. In 
this case the cavity is in the vacuum state with zero photons (n = 0), while the 
atom is prepared in its excited state |e). The probabilities for the naked basis states 
{|e, 0), \g, 1)} are given by Eqs. (|57|) with n = resulting in 

P e ,o(t) = cos 2 (^t) = \ [1 + cos(^ot)] 

(52) 

P gA {t) = sin 2 (^t) = i [1 - cos(^ t)] 
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where i?o = 2g is the vacuum Rabi frequency, that is the transition frequency for the 
dressed states | ± 0) = (|e, 0) ± \g, l))/\/2 in the subspace 5(0), with energies 

E±o = h{u±^) = H(u±g). (53) 

The energy separation between the two dressed levels is (AE)q = 2hg and this par- 
ticular situation is known as vacuum Rabi oscillations or 1-photon oscillations. This 
means that the spontaneous emission of a photon by the atom inside a cavity becomes 
a reversible process, with an oscillatory behavior in time, which is in hard contrast 
with the irreversible spontaneous emission in free space, with exponential decay in 
time. This is due to the radically different energy spectral distribution of the electro- 
magnetic field in the two cases: in the first there is only one field mode in the cavity 
accessible by the emitted photon, in the second there is a continuum of available 
modes. On the other hand, the spontaneous emission in a cavity can be suppressed 
if there is not any field mode resonant with the atomic transition. Modifying the 
field spectral distribution the spontaneous emission can be enhanced or suppressed^ 
(cavity QED). 

A simple explanation of the vacuum Rabi oscillatons can be achieved by thinking 
the atom (weakly excited since the field mode has a maximum of one photon) as an 
oscillating dipolc, in particular a quantized harmonic oscillator. The interacting atom 
and mode can be thought as two weakly coupled oscillators at the same frequency, 
described by the Hamiltonian 

n JC ->U = fiw{tfb + at a) + Tig(a ! b + atf) (54) 

where the commutation relations holds for the two subsystems [a, a'] = 1 (mode) and 
[w , b] = 1 (atom). This Hamiltonian can be diagonalized introducing the two normal 
modes operators 

obtaining the new one 

% = H(uj + g)c i c + h(uj - g)dU. (56) 

This is the Hamiltonian describing two independent harmonic oscillators (normal 
modes) of different frequencies u> ± g. 
The eigenstates of H are 

H\n c ,n d ) = [H(oj + g)n c + h(u) - g)n d ]\n c ,n d ) (57) 

which become for the lowest states 

H\l c ,Q d ) = h(Lu + g)\l c ,0 d ) ; H\0 C , 1 d ) = h(u> - g)\0 c , l d ). (58) 

In terms of the original mode operators a and b, the eigenstates can be compared to 
the dressed states for the atom and the cavity field mode: 

'|lc, d ) = ^ (\l a ,0 b ) + \0 a , 1 6 )) f | + 0) = X (\e, 0) + \g, 1)) 

V2 ^ { V2 (59) 

|0« U) = ^ (|l o ,0i) - |0 OJ 1 6 » ( | - 0) = ^ (|e, 0) - \g, 1)) 



3 The spontaneous emission can be also suppressed in dielectric media with periodically 
time changing refractive index (photonic crystals), where there are energy bands "forbidden" 
for photons, just like the case of electronic bands in crystals. 
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Hence, the lowest stationary states of two harmonic oscillators correspond to the 
energies K(w ± g), exactly as the first two dressed states | ± 0) of the lowest manifold 
~(0) of the JC model. 

Experimental confirmations of Rabi oscilations were given by the group of Kimble 
[21) . in the optical domain, and by the group of Haroche [5], in th microwave domain. 
The atom-field system in the vacuum state can be weakly excited by a probe field 
tunable in frequency. The outcome is that the excitation is not resonant with the field 
mode (or the atomic) frequency ui, but, as predicted, with the two frequencies of the 
normal modes (c and d) of the coupled system: w ± g. 

8.1 An application: swapping of the excitation 

Suppose that a two-level atom, prepared in the superposition state c e |e) + c g \g), is 
injected in a cavity in the vacuum state |0), giving as initial state for the whole system 

\*(0)) = (c e \e)+c a \g))\0). (60) 

Using Eqs. (|4"3")h one derives the following simple time evolutions for the states |e,0) 
and Is-, 1): 

|e,0) ->■ cos(^) |e,0) -isin(^t) \g, 1) 

\ g> i) C os (%t) \g, 1) - isin (Ait) |e, 0) (61) 

|S, 0) -> \g, 0) (uncoupled) 

If the system evolves for a time t such that f2ot = ir the above system of equations 
give the results 

'|e,0)^-%,l) 

< Iff, 1> -»> — i|e, 0> (62) 

>o>->|ff,o) 

and the whole system starting from the state (|60p is in the state 

|^(t)) = | 5 )(-ic e |l)+c a |0)) (63) 

from which it is evident that the initial superposition state of the atom has been 
swapped with (or transferred to) the state of the cavity mode (except for a local 
phase factor). 

Since the process is reversible, one may imagine that after the time t the atom goes 
outside the cavity and another one is injected in the ground state | (7)2 to interact 
with the cavity mode (now the state of the whole system is given by Eq. (|6"3")l ). After 
a second time evolution with figt — ir the outcoming state of the system is 

|ff)a(-ice|l) + cg\0)) -> (-c e |e) 2 + c g \g) 2 )\0) (64) 

which means that the photonic state of the cavity now has been swapped with the 
second atom, recovering a state similar to the initial one. 

9 Collapses and revivals 

Going back to the general case of n 7^ (always for resonance conditions), it has been 
shown that if the field mode is prepared in a Fock state \n) and the atom injected in 
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the cavity in the state |e), the probability to find the atom again in its excited state 
after a time t > is given by the first of Eqs. ((37)) . 

If the field state is prepared in a generic superposition of Fock states pf{0) = 
2^°=o -Pnl n )( n |j the probability to find the atom in the excited state is 

oo * 

p e (t) = J2pn p e,n(t) = 2 

n=0 

which is an average over all the possible (infinite) Fock states \n) weighed by the 
probabilities P n . 

If, in particular, the initial state of the cavity is a thermal state with P n T ^ = 




Fig. 5. Collapse of Rabi oscillations for the cavity field initially in a thermal state with 
(iV) = 10 and the atom in the excited state |e). 



oo 

^ Pn COS 

n=0 



(65) 



(i+(jv))"+ 1 ' the infinite sum of sinusoidal oscillations, each weighed by P n , brings 
the state to a rapid collapse of Rabi oscillations, until P e (t) \ with some small 
fluctuations (see Fig. [3]). In this way the probabilty to observe the atom in its excited 
or ground state is equally distributed, exhibiting a "classical" effect. 
Different is the case of an initial coherent state |a) of the cavity field. If the co- 
herent state amplitude is too high, the field is "classical" and, as for the Fock state 
\n), there is ideally a Rabi oscillation at a mean frequency J?/„\ ~ 2g\a\ — 2g^ (n) a 

(semi-classical thoery, see Ramsey cavities). For a coherent state P n C ^ = e ~H 2 1°^ 
and the atomic excitation probability is 



Pe(t) = 



1 



1 



x >a\ 2n 



cos(fi n t) 



n=0 



(66) 



First of all it is possible to demonstrate that for \a\ 2 3> 1 the probability for the 
atomic excited state is 



PS) 



1 + cos(25|a|i)e"' a ^" 



(67) 



where the oscillations occur at the mean Rabi frequency £2i n \ — 2g^J (n) a = 2g\a\ 
with a Gaussian envelope of width t c ~ ^ , known also as coherence time (or collapse 
time). This result can be interpreted in terms of a destructive interference of Rabi 
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Fig. 6. Collapses and revivals of Rabi oscillations for the cavity field initially in a coherent 
state | a) with |a| 2 = 25 and the atom in the excited state |e). In evidence there are three 



of the revival times gt 
central frequency i?(„) . 



j2Tr\a\ and in the inset it is shown the Rabi oscillation at the 



oscillations at different frequencies, where only a narrow range of Rabi frequencies 
Afl « g around Qi n \ plays a role in the dynamics (in contrast with the thermal case). 
The mechanism of collapse is always interpretable as a "classical" effect. 
What is a purely quantum feature is the revival effect of Rabi oscillations, which 
occurs with a period of time 

„ M 

t r ~ Z7T > t c . 

9 

The revival time is the period of time for which two Rabi oscillations at two close 
frequencies (and necessarily close to the central frequency) provide a constructive 
interference. This happens when the two oscillations get in phase and an extimate of 
the revival time is the following 

t V) = M = ~ (68) 

r %-%)_! 2g(\a\-y/\aF=T) 9 

for \a\ 2 3> 1, but not too high. The revival effect, shown in Fig. [5] is a phenomenon of 
quantum recurrence due to the "granular" nature of the quantized electromagnetic 
field and it allows the oscillations at different frequencies to be periodically in phase. 
For long times the effect gets more and more ruined since the revivals are never 
complete. 



10 The Haroche experiment 

In 1996 the group of Haroche, from the Ecole Normale Superieure of Paris, published 
a paper [5] named "Quantum Rabi Oscillation: A Direct Test of Field Quantization in 
a Cavity" . The realization of the experiment required a resonant interaction between 
a two-level atom and one cavity mode, a strong coupling between the two subsystems 
with an interaction time shorter than the atomic decay, a high-Q cavity, the control 
of paramaters like atomic excitation and speed, and cavity temperature. A supercon- 
ductive microwave cavity was chosen, through which the atoms, selected in velocity 
and excited at Rydberg levels (main quantum number n ^> 1), passed one at a time, 
with one Rydberg transition resonant with the cavity mode. 
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Fig. 7. Results of the Haroche experiment [5], which put in evidence collapses and re- 
vivals of Rabi oscillations, depending on the mean photon number present in the cavity: (A) 
{n) a ~ 0, (B) (n) a ~ 0.40, (C) (n) a ~ 0.85 and (D) (n) a ~ 1.77. (a), (b), (c), (d) are the 
corresponding Fourier transforms, whereas (a), (/?), (7), (<5) are the corresponding photon 
number distributions. 



In particular the group of researchers used Rubidium atoms excited in circular Ryd- 
berg states (|e) -O- n e = 51 and \g) o n g = 50), with a transition frequency between 
the two states v eg ~ 51GHz. The average lifetime of the atomic excited level was 
77 ~ 32ms, whereas the interaction time between the atom and the cavity mode 
was ti n t — 10ms. Using Rydberg atoms ensures that the dipole coupling (d oc n 2 ) is 
very strong, more precisely the coupling constant was g ~ 1.6 x 10 5 Hz. The supercon- 
ductive microwave cavity {L ~ 3cm) supported one mode resonant with di Rydberg 
transition v = v eg . The cavity was cooled at a teperature of T = 0.8K, corresponding 
to n ~ 0.05 thermal photons. The quality factor of the cavity was Q — ~ 108, 
which means that the average lifetime of cavity photons was ~ 220/xs. The values 
of these parameters allowed a strong coupling regime, for which the coupling frequency 
was much greater than all the decay rates g 3> {2k, 7,4} and a coherent JC dynamics 
may occur. Moreover a source of microwave radiation could inject into the cavity a 
small coherent field (from zero to few photons) with a high controlled energy. Detec- 
tors based on field ionization techniques allowed to count the outcoming atoms in the 
excited or in the ground state. The atoms were prepared, injected and measured one 
at a time, separated by a reset time (= 2.5ms). It was possible to obtain the occu- 
pation probabilities P e and P g in function of the interaction time t- mt , which could 
be changed varying the atomic velocities. After reaching a sufficient statistics, the 
ensemble of measurements was repeated for different interaction time values. These 
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measurements referred to the probability 

1 - y^P n cos(i?„t) 

with t = ti nt (Fig.[7R-D), and also to the Fourier transform of P g (t) (Fig.[7^-d) and to 
photon statistics (Fig. [7]a:-(> j3- All these measured quantities were repeated for differ- 
ent coherent field amplitudes of the cavity mode state. In the first case (n) a = 0, the 
vacuum Rabi oscillations were observed, for the first time in the time domain, at the 
frequency v$ = S 3 - ~ 47kHz. For a greater coherent field amplitude inside the cavity 
(n) a = 0.4,0.85, 1.77 the signal P g (t) lost progressively its sinusoidal behavior, with 
an evidence of the presence of collapses and revivals in Fig. [7p-D. From the Fourier 
transforms plots (a-d) there is evidence of subsequent Rabi frequencies involved in 
the dynamics, as the coherent field amplitude increases, in particular the frequencies 
v$yj\2), voy/~[3), 2vq as predicted by the theory v n = 5? = v^^/n + 1. This is the 
witness of the granular nature of the quantized electromagnetic field inside the cav- 
ity. The weights of the Fourier components yielded the photon number distributions 
in the field (a-S). For no injected field the pionts are in good agreement with the 
theoretical exponential decay curve of the thermal radiation. When a coeherent field 
is present in the cavity, the photon number distribution follows a Poisson statistics, 
signature of a coherent radiation. 



P g (t) = l-P e (t) = - 



11 Generation of two atoms entangled states 

Another important experiment from the group of Haroche [B], published in 1997 

with a paper named "Generation of Einstein- Podolsky- Rosen Pairs of Atoms" , based 

essentially on the same apparutus as the previous one, exploited the resonant JC 

interaction to generate entangled states of two atoms via a cavity mode. 

At time t — a first circular Rydberg two-level atom was prepared in the excited 

state |e)i and injected in a superconductive microwave cavity in the vacuum state 

|0>: 

|f(0)) = |e)i®|0) = |ei,0). 

After a certain interaction time t\ such that fl{)t\ — -|, where fl§ = 2g is the vacuum 
Rabi frequency, the first atom interacted with the cavity field mode. Remembering 
Eqs. (l6"TT) in Section IQ1 and using real coefficients, one has 

|e, 0) ->■ cos (fit) |e, 0) + sin (ty-t) \g, 1) 

(69) 

\g, 1) -> cos \g, 1) + sin (^t) \e, 0) 
which brings to evolved state 

|e 1 ,0)->|!?(t 1 )) = i(|e 1 ,0) + | Sls l)). (70) 

After a time interval T a second atom, prepared now in the ground state |<72), was 
injected in the cavity and the state is 

\W{tt +T)) = \g 2 ) ® (|ei,0) + \gi,l)) = -j= • 



4 Reprinted figure with permission from M. Brune et ai, Phys. Rev. 
Lett. 76, 1800 (1996). Copyright (1996) by the American Physical Society. 



http://link.aps.org.pros.lib.unimi.it/abstract/PRL/v76/pl800 
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The second atom interacted with the cavity mode for a time t 2 — 2ti, such that 
n t 2 = 7r and, since the state \g 2 , 0) is uncoupled, the state of the system was 

+ T + t 2 )) = \^,0) + \ 9l ,e 2 ,0) = w _ Uia2 ^ |q)f 

where \&~}aiA2 = l ei ' g2 ^ gl ' e2 ^ is one of the maximally entangled states of the Bell 
basis. The cavity field was in the vacuum state and it "catalysed" the process of en- 
tangling two independent atoms through the common and subsequent JC interaction. 
After a time t > t\ + T + t 2 the energy of the two atoms was measured or, in order 
to confirm the generation of the entangled state, interference effects measurements 
(Ramsey fringes) or correlation measurements (Bell signal) were performed. 



12 Conclusions 

This tutorial intends to provide the basic elements for the comprehension of the 
Jaynes-Cummings model, in order to give a pedagogical approach to the study of 
matter-radiation interaction with the tools of quantum mechanics. Along with this 
approach, the Hamiltonian of the interacting system is derived from general state- 
ments and its eigenvalues and eigenvectors are calculated from diagonalization. The 
power and simplicity of the JC model stem from the fact that the dynamics occurs 
within a bidimensional subspace spanned by the dressed basis states, which are re- 
lated to the naked basis states by a simple rotation matrix. The description of the 
dynamics by means of these two basis explains in an elegant way the Rabi oscillations 
phenomenon, that is the coherent exchange of energy between a two-level system and 
a resonant electromagnetic field mode. A particular attention is devoted to the vac- 
uum Rabi oscillations, where the initial state of the field mode contains zero photons. 
The formal description of the evolution operator is presented in order to obtain the 
JC dynamics for every initial state of the system, including the time evolution of the 
subsystems. The real quantum feature of the JC model is given by the revival effect, 
occurring with an initial coherent state of the field mode, which can be explained 
only if the electromagnetic field is quantized. 

cQED is a natural framework for the implementation of the JC model, where single 
two-level Rydberg atoms can pass through or be trapped in microwave single mode 
cavities. Two key-experiments are described in the last sections, where Rabi oscilla- 
tions, together with their collapse and revival effects, and the swapping application 
for retrieving maximally entangled states, are the results of spectacular observations 
and measurements. 
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